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This  research  project  is  concerned  with  the  development  of  novel  nano-engineered  materials 
whose  thermal  properties,  such  as  thermal  conductivity,  can  change  in  a  predetermined  manner 
due  to  internal  and  external  triggers,  such  as  temperature  and  mechanical  stresses.  The  capability 
to  control  actively  or  passively  the  thermal  properties  of  materials  will  allow  for  adapting 
dynamically  the  thermal  transport  paths  to  time  varying  operating  conditions  and  enable  radically 
improved  thermal  management  systems,  allowing  for  larger  power  densities  in  electronic 
components  and  larger  operating  temperatures  of  load  barring  structures  as  well  as  reducing 
thermally  induced  damage.  The  mechanisms  that  can  dynamically  change  thermal  properties  of 
materials  are  dominated  by  submicron  and  nanoscale  phenomena.  In  order  to  design 
systematically  these  materials,  their  submicron/nanoscale  behavior  needs  to  be  modeled  and 
predicted  with  sound  accuracy  but  time  and  cost  maneuverable. 

With  the  AFOSR  support,  the  Pis  have  developed  a  novel  computational  framework  for  the 
analysis  and  design  of  nano-structured  materials  with  tunable  thermal  properties.  The 
accomplishments  of  this  3-year  project  including: 

1) .  A  new  multiscale  modeling  approach  based  on  the  Boltzmann  transport  equation  (kinetic 
theory)  for  predicting  phonon  transport  has  been  developed.  This  computational  tool  predicts 
adequately  the  thermal  performance  of  the  nanostructures  balancing  accuracy  and  turn-around 
time,  and  has  been  verified  through  comparison  with  standard  Monte-Carlo  simulation  methods 
and  experimental  studies. 

2) .  Frequency-dependent  phonon  transmission  across  interfaces  of  dissimilar  materials  has  been 
simulated  using  the  nonequilibium  Green’s  function  (NEGF)  approach.  Size-dependent  phonon 
transmission  has  been  studied  for  material  interfaces  with  transverse  or  longitudinal 
confinements.  This  size-dependent  phonon  transmission  is  the  key  input  parameter  for  the 
Boltzmann-equation  based  multiscale  modeling  tool  for  nanostructured  materials. 

3) .  We  explored  the  dynamical  tuning  of  Mthermal  conductivity  of  materials  using  mechanical 
stress/strain.  We  have  shown  that  strain  effects  on  the  lattice  thermal  conductivity  of  inorganic 
nanostructures  and  bulk  polymer  materials  can  be  significant  through  molecular  dynamics  (MD) 
simulations. 

4) .  An  optical  pump  and  soft  X-ray  probe  measurement  system  has  been  developed  to  measure 
nanoscale  thermal  transport.  Our  work  made  the  first  observation  and  quantitative  measurement 
of  non-Fourier  heat  dissipation  from  a  nanoscale  hot  spot  to  the  surrounding,  which  was  recently 
published  in  Nature  Materials  in  2010. 

We  have  published  2  book  chapters  and  8  papers  in  high  impact  journals  in  2008-2010  (1  in 
Nature  Materials,  4  in  Phys.  Rev.  B,  2  in  J.  Appl.  Phys.,  1  in  Int.  J.  Num.  Meth.  Engr.).  In 
addition,  we  currently  have  3  journal  papers  accepted  and  5  journal  papers  under  review  as  of 


February  2011.  In  addition,  we  have  given  about  30  invited  and  contributed  presentations 
acknowledging  the  AFOSR  support  for  this  project. 


In  the  following,  we  highlight  the  afore-mentioned  achievements: 


1).  Multiscale  Thermal  Transport  Modeling  based  on  Boltzmann  Transport  Equation 


Although  the  Boltzmann  transport  equation-based  kinetic  formulation  is  an  attractive  modeling 
and  analysis  approach  for  a  widely  applicable  computational  framework,  solving  the  7- 
dimension  (3  in  space,  2  in  direction,  1  in  frequency  and  1  in  time)  integro-differential  type 
Boltzmann  equation  is  a  prohibitively  expensive  and  time-consuming  task.  Existing  studies  for 
phonon  transport  in  nanoscale  composites  [1,  2,  3,  4]  have  several  common  assumptions  that 
may  not  be  appropriate:  the  shape  of  nanowires  or  nanoparticles  in  existing  studies  is  regular, 
such  as  circular,  rectangular  or  cubic,  the  size  of  the  nanoparticles  and  nanowires  is  typically 
assumed  to  be  uniform,  and  the  distribution  of  these  nano-inclusions  is  normally  assumed  to  be 
periodic.  Careful  examination  of  the  SEM  images  [Figure  1(a)]  clearly  shows  that  neither 
assumption  is  true  for  realistic  nanocomposites.  Extending  either  the  discrete  ordinate  method  or 
Monte  Carlo  simulation  to  more  general  nanostructures  is  not  trivial.  To  analyze  heat  conduction 
at  submicron  scales  of  geometrically  complex  nano-structured  materials,  we  have  developed  an 
extended  finite  element  method  (XFEM)  is  to  solve  the  phonon  Boltzmann  transport  equation  as 
shown  in  Fig.  1(b)  and  (c)  [5].  The  geometry  of  the  material  interfaces  is  described  by  a  level-set 
approach.  Discontinuities  of  the  phonon  distribution  across  material  interfaces  are  captured  via 
enriched  shape  functions.  To  enforce  interface  scattering  conditions  and  boundary  conditions,  a 
stabilized  Lagrange  multiplier  method  is  used.  The  utility  of  the  XFEM  approach  is 
demonstrated  through  the  thermal  analysis  of  experimentally  characterized  material  samples  [See 
Figure  1(c)],  The  XFEM-BTE  could  potentially  be  extended  to  study  thermoelectric  and 
electrical  properties  of  nanocomposites  when  electron  BTE  is  solved  along  with  the  Poisson 
equation. 


(a) 


Figure  1.  extended  finite  element  method  for  phonon  Boltzmann  transport  equation  (XFEM-BTE)  has  been 
developed  to  study  thermal  conductivity  of  nano-enabled  bulk  materials  with  complex  structures:  (a).  TEM  image  of 
the  ball  milled  SiSoGe2o  nano  composites  [6].  (b).  Geometry  described  by  level-set  approach,  (c).  Temperature 
contours  simulated  with  XFEM-BTE. 


In  general,  a  large  number  of  acoustic  and  optical  phonon  distributions  need  to  be  simulated  to 
accurate  capture  the  dispersion  characteristics  of  the  individual  material  phases.  In  this  project,  a 
simpler  and  numerically  less  costly  approach  was  pursued  based  on  the  single-distribution  quasi- 


static  gray  phonon  model.  Approximating  the  collision  operator  by  a  BGK  approach,  the 
governing  equations  of  the  gray  phonon  Boltzmann  transport  equation  (GPBTE)  are: 

e°  -p 

s- Ve  = -  (x,p)eQxS  (1) 

A 

where  the  vector  s  denotes  the  propagation  direction  of  the  phonon;  e  is  the  phonon  distribution, 
e°  the  equilibrium  distribution  and  A  the  phonon  mean  free  path  length.  Note  the  phonon 
distribution  is  defined  as  function  of  the  spatial  domain  Q  and  the  velocity  space  S  with  x 
denoting  the  spatial  and  p  the  velocities  coordinates. 

Assuming  an  infinitely  periodic  material,  the  GPBTE,  eqn.  (1),  is  solved  in  a  unit  cell  imposing 
proper  periodic  and  symmetry  boundary  conditions.  To  capture  phonon  scattering  at  material 
interfaces,  a  diffusive  mismatch  model  (DMM)  has  been  adopted.  More  complex  interface 
models  can  be  easily  integrated  into  the  numerical  framework  developed  in  this  project. 

To  represent  the  geometry  of  the  material  interfaces,  a  level-set  presentation  of  the  material 
layout  is  used.  The  level-set  representation  is  reconstructed  from  a  two-dimensional  image  of  a 
material  sample  using  standard  image  processing  techniques.  The  material  interface  is  defined  by 
the  zero  level-set.  To  avoid  the  need  for  body- fitted  meshes,  an  extended  finite  element 
approach  (XFEM)  has  been  adopted.  This  approach  allows  capturing  the  discontinuities  at  the 
material  interfaces  on  structured  meshes.  A  discontinuous  Galerkin  method  has  been  identified 
as  the  method  of  choice  for  discretizing  the  velocity  space. 

XFEM  has  been  originally  developed  for  modeling  crack  propagation  in  solids  which  leads  to  a 
discontinuity  in  the  displacements  at  the  crack  boundary.  More  recently  XFEM  has  been  applied 
to  a  broad  range  of  problem  involving  discontinuous  fields,  such  as  phase  change  problems  and 
multi-phase  flows.  In  this  project  XFEM  was  adopted  for  solving  the  GPBTE  with 
discontinuities  in  the  phonon  distributions  along  the  material  interfaces.  In  addition  to  the 
standard  Galerkin  approximation  of  the  phonon  distribution  a  discontinuous  enrichment  function 
is  added: 

e(x,p)  =  YjNiP  (x)  [e,00  +  ^(x)  d,(p)]  (2) 

i 

where  the  ^(x)is  the  enrichment  function  and  d;.(p)  the  associated  unknown.  To  capture  the 

discontinuity  at  the  phonon  interface,  a  Heaviside  function  is  used  to  enrich  the  phonon 
distribution,  i.e. 

f +1  for  xeQ, 

^(x)  =  tf(x)  =  J  (3) 

[0  for  x  e  Q7/ 

Note  the  enrichment  is  only  active  in  elements  intersected  by  a  material  interface,  i.e.  the 
unknown  n,  (p)  is  set  to  zero  in  all  elements  that  are  not  intersected  or  adjacent  to  an  intersected 
element. 

The  approximation  (2)  is  used  to  discretize  in  space  the  weak  form  of  the  governing  equations 
(1).  As  the  problem  is  convection  dominant,  the  discretized  problem  suffers  from  numerical 
instabilities  and  needs  to  be  stabilized.  In  this  project  a  Streamline-Upwind  Petrov-Galerkin 
(SUPG)  stabilization  approach  is  used.  In  the  context  of  GPBTE  the  numerical  diffusion 
introduced  by  SUPG  can  be  well  controlled  and  leads  to  minor  additional  computational  cost. 

To  enforce  the  interface  conditions  of  the  DMM  a  stabilized  Lagrange  multiplier  formulation 
was  developed.  To  satisfy  the  inf-sup  stability  condition  for  mixed  finite  element  formulations,  it 
is  usually  sufficient  to  select  the  approximation  space  for  the  Lagrange  multipliers  to  be  one 


polynomial  order  lower  than  the  space  for  the  physical  field,  i.e.  the  phonon  distribution.  To 
overcome  the  additional  complexity  due  the  stabilized  Lagrange  multiplier  formulation,  a  much 
simpler  phonon-flux  interface  condition  was  recently  developed.  This  formulation  is  derived  via 
integration  by  parts  the  weak  form  of  the  governing  equations  (1).  For  example  for  phase  I,  this 
leads  to: 


s-We  +—(eI  -e/0) 
AV  ’ 


dQ  =  ^5e s-n  e1  dT -  J SVe ■  s  e1  dQ.  +J 8e ~-(g/  - )  dQ. 


(4) 

The  boundary  term  allows  satisfying  the  DMM  conditions  without  Lagrange  multipliers  by 
substituting  the  incoming  fluxes  with  the  DMM  condition  (2)  as  follows: 

j  de  s  -n  e1  dT  =  J  Se  s  ■  n  e1  dT  +  J  Se s-ne1  dT 


sn>0 


sn<0 


(5) 


with 
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In  this  project,  the  XFEM  formulation  of  the  GPBTE  was  studied  extensively  for  a  discrete 
ordinate  discretization  of  the  velocity  space.  Reference  results  were  taken  from  the  square- 
inclusion  example  of  Chen  and  Yang  [7]  (see  figure  2).  Additional  problems  using  real  samples 
of  nano-composites  were  studied.  Numerical  studies  suggested  that  the  XFEM  formulation 
features  the  same  accuracy  and  stability  as  finite  element  and  finite  volume  formulation  requiring 
body-fitted  meshes.  Comparison  between  XFEM  predictions  and  experimentally  measured  heat 
conductivities  from  Joshi  et  al.  2008  showed  good  agreement  [see  Figure  1], 


o 

Q 

o 

~o 

O 

’s_ 

CD 

Q. 


23.5  nm 


I  23.5  nm 


o 

Q 

u 

“O 

O 

‘i_ 

cu 

Q. 


0.8 1 
0.6 
0.4 
0.2 
0.0 
-0.2 


-0.4 


-0.6 


-0.8 

0. 


*  y/L=0.50 
□  y/L=0.70 
o  y/ 1=0. 85 


symmetry  be  Temperature  contours 

Figure  2:  XFEM  Analysis  of  synthetic  nano-comoosite  wi. 
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However,  a  significant  limitation  of  our  current  XFEM-BTE  is  that  we  used  the  acoustic 
mismatch  model  (AMM)  [8],  and  diffuse  mismatch  model  (DMM)  [8]  for  phonon  transmission 
across  material  interfaces,  which  are  standard  practice  today  by  various  groups  but  are  indeed 
incorrect  for  realistic  material  interfaces. 


2  Phonon  Transmission  across  Coherent  Interfaces 

Realizing  that  carrier  transport  across  interfaces  is  the  key  to  the  superior  properties  of  nano- 
enabled  bulk  materials  and  the  currently  available  AMM  and  DMM  model  could  not  be  valid  in 
realistic  materials,  we  have  investigated  a  variety  of  modeling  approaches  for  phonon 


transmission  across  interfaces.  Molecular  dynamics  simulation  has  been  used  extensively  to 
study  the  interface  thermal  resistances  of  various  materials.  However,  the  results  could  only  be 
used  to  qualitatively  understand  the  impact  of  interfaces,  but  not  useful  for  qualitative  design, 
since  it  cannot  be  scaled  up  or  integrated  into  a  multiscale  simulation  tools. 

The  multiscale  BTE  solvers  described  in  session  1)  essentially  need  frequency-dependent 
phonon  transmission  or  energy-dependent  electron  transmission  across  interfaces  as  input 
parameters.  For  phonon  transmission  across  interfaces,  atomistic  methods  such  as  phonon  wave- 
packet  method  [9],  linear  lattice  dynamics  approach  [10]  and  nonequilibrium  Green’s  function 
(NEGF)  approach  [11-20]  has  recently  been  developed.  Among  these  atomistic  methods,  NEGF 
is  the  most  efficient  one  for  calculating  phonon  transmission  across  interfaces  for  considering 
detailed  atomic  micro-structures.  The  NEGF  approach 
calculates  phonon  transmission  by  solving  the  lattice 
dynamics  equation  under  harmonic  approximation 
using  the  Green’s  function  method. 


Under  harmonic  approximation,  the  dynamical 
equation  of  a  lattice  system  can  be  solved  with  the 
Green’s  function  method, 

(cq2\-H)G  =  1,  (6), 

where  H  is  the  dynamical  matrix  and  G  is  the  Green’s 
function.  The  Green’s  function  represents  the  response 
of  the  system  to  an  infinitesimal  perturbation.  For  a 
general  interface  system  which  consists  reservoirs  1 
and  2,  and  an  interface  region  as  shown  in  Figure  3,  the 
total  phonon  transmission  across  the  interface  region 
can  be  formulated  as, 

’E(oj)  =  Trace[r]GDDr2G+DD],  (7), 

where  GD  D  is  the  Green’s  function  of  the  interface  region  solved  from  equation  (6)  and  T1>2  are 

related  to  the  self-energy  matrices  due  to  the  coupling  between  the  reservoirs  and  the  interface 
region,  which  can  be  calculated  from  dynamic  matrix  Hx  2  and  the  interaction  r,  2 . 

The  transmission  function  at  each  particular  phonon  frequency  can  be  written  as, 

t(a>)  =  E(©)  /  SO)Reservoirl ,  (8), 

where  S(&>)Reservoirl  is  the  total  phonon  transmission  in  pure  materials  1. 

The  above  NEGF  method  is  a  general  method  and  can  be  applied  to  any  atomistic  system  under 
harmonic  approximation.  However  the  application  of  this  method  can  be  greatly  limited  due  to 
the  challenge  of  solving  equation  (6)  if  the  system  size  is  too  large.  To  make  the  matrix  system 
solvable,  the  studies  so  far  using  NEGF  usually  assume  that  the  atomistic  system  is  lattice 
matched  (coherent  interface  with  same  lattice  constants).  A  number  of  systems  have  been  studied 
with  NEGF  approach,  such  as  low  dimensional  molecular  chains,  coated  nanowires,  amide- 
linked  nanotubes  and  lattice  matched  Si/Ge,  etc. 

Figure  4  shows  our  recent  studies  [X.B.  Li,  et.  al.,  Submitted  to  PRB],  Figure  4(a)  shows  that 
phonon  transmission  across  Si-Ge  interface  (assuming  same  lattice  constant  of  the  two  materials) 
calculated  using  NEGF  agrees  well  with  the  one  obtained  from  lattice  dynamics  [10].  Figure  4(b) 
shows  phonon  transmission  across  Si/Ge  like  materials  interface  by  varying  the  mass  in  Ge-like 
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Figure  3.  A  general  interface  system  for  phonon 
transmission  calculation  using  NEGF  approach. 
The  system  is  generally  divided  into  three  parts: 
reservoir  1,  interface  region  and  reservoir  2.  H  is 
the  dynamical  matrix  of  each  part  and  i 
represents  the  interactions  between  the  reservoirs 
and  the  interface  region. 


material  while  keeping  the  Ge  force  filed.  With  varying  mass  ratio  phonon  dispersion  mismatch 
can  be  tuned  and  significant  change  can  be  obtained  for  phonon  transmission.  Figure  4(c)  shows 
the  thermal  conductance  across  the  Si/Ge  like  materials  interface  and  an  optimal  mass  ratio  at 
1.26  exists  for  achieving  the  maximum  thermal  conductance. 
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Figure  4.  (a)  Phonon  transmission  across  lattice  matched  Si-Ge  interface  using  NEGF.  The  results  agree  well  with 
lattice  dynamics  calculations,  (b)  phonon  transmission  across  S(MI)/Ge  like  (M3)  interfaces  with  varying  mass 

ratio:  CC  =  Af  3  /  M] .  (c)  Thermal  conductance  across  the  interface  as  a  function  of  the  mass  ratio. 


However,  the  current  NEGF  approach  has  assumed  lattice-matched  systems  and  harmonic 
approximations,  which  greatly  limits  the  power  of  NEGF  approach.  For  example,  strain  effects 
on  phonon  transmission  has  been  studied  for  the  lattice  matched  Si/Ge  interface  [14]  and  the 
results  only  show  very  small  difference  from  the  unstrained  interface,  which  does  not  agree  with 
experimental  results  [21].  This  suggests  the  importance  of  considering  the  realistic  material 
interfaces  (with  lattice  mismatch).  The  anharmonic  effects  have  also  been  addressed  in  several 
studies;  however  the  method  has  only  been  applied  to  simple  one-dimensional  chains  due  to  the 
computational  complexity.  Further  developments  need  to  be  focused  on  addressing  these 
challenges  in  realistic  material  interfaces. 

3.  Tuning  Thermal  Conductivity  of  Materials  using  mechanical  Stress/Strain. 

Different  methods  can  be  used  for  tune  thermal  conductivity  of  materials.  Among  these  methods, 
mechanical  strain  provides  an  efficient  way.  We  have  performed  molecular  dynamics  (MD) 
simulation  to  systematically  study  the  strain  effects  on  the  lattice  thermal  conductivity  of 
inorganic  materials  (mainly  nanostructures:  Si  nanowire  and  thin  film,  single-walled  carbon 
naotube  (SWCNT)  and  single  layer  graphene)  and  bulk  polymer  materials  [X.B.  Li,  et.  al.,  PRB 
(2010),  J.  Liu,  et.  al.,  PRB  (2010)].  Results  show  that  thermal  conductivity  of  the  strained  silicon 
nanowires  and  thin  films  decreases  continuously  when  the  strain  changes  from  compressive  to 
tensile.  However,  the  thermal  conductivity  has  a  peak  value  under  compressive  strain  for 
SWCNT  and  at  zero  strain  for  single  layer  grapheme  as  shown  in  Fig.  5(a)  for  results  on 
SWCNT.  In  contrast,  thermal  conductivity  of  polymer  materials  increases  with  increasing  tensile 
strain  as  shown  in  Fig.  5(b).  The  underlying  mechanisms  are  analyzed  in  this  paper  for  both 
types  of  materials.  We  found  that  the  thermal  conductivity  of  inorganic  materials  can  be  related 
to  the  phonon  dispersion  curve  change  and  structural  buckling  under  strain  and  for  polymer 
chains  thermal  conductivity  directly  connects  to  the  orientations  of  the  chains  instead  of  strain 
rates.  This  thermal  conductivity  dependence  with  strain  can  guide  us  to  tune  the  thermal 
conductivity  for  materials  in  applications. 


Figure  5.  (a).  Strain  effects  on  the  lattice  thermal  conductivity  of  armchair,  zigzag  and  chiral  SWCNTs.  (b).  Strain 
effects  on  bulk  polymer  materials.  Different  strain  effects  are  observed  for  these  two  kinds  of  materials. 


4.  Characterizing  Nanoscale  Thermal  Transport  using  Soft  X-Ray  Probes 

Femtosecond  laser  is  a  unique  tool  to  study  a  number  of  ultrafast  relaxation  processes  and 
nanoscale  transport  phenomena.  Fourier  theory  of  heat  conduction  considers  heat  transport  as  a 
diffusive  process  where  energy  flow  is  driven  by  a  temperature  gradient.  However,  at  length  scales 
smaller  than  the  mean  free  path  for  the  energy  carriers,  heat  flow  becomes  ballistic  -  driven  by 
direct  point-to-point  transport  of  energy  quanta.  Although  past  experiments  have  seen  ballistic 
effects  in  layered  thin  films  or  nanowires,  non-Fourier  heat  dissipation  from  a  nanoscale  hot  spot  to 
the  surrounding  had  not  been  experimentally  observed  and  was  a  subject  of  some  uncertainty 
although  it  has  significant  relevance  in  heat  dissipation  of  electronic  transistors.  Our  work  made 
the  first  observation  and  quantitative  measurement  of  this  geometry  using  the  excitation  and 
detection  scheme  shown  in  Fig.  6  in  exciting  recent  work  recently  published  in  Nature  Materials  in 
2010  [M.  Siemens,  et.  al,  Nature  Materials  (2010)].  We  use  a  beam  of  ultrafast  coherent  soft  x-rays, 
at  a  wavelength  of  29  nm,  to  directly  observe  the  cooling  dynamics  of  a  nanoscale  heat  source  into 
its  bulk  surroundings.  By  interferometrically  monitoring  displacement  in  a  heated  nanostructure 
using  diffraction  of  soft  x-ray  light,  we  can  detect  dynamic  temperature  changes  more  sensitively 
than  would  be  possible  using  either  ultrashort  optical  pulse  probing  or  steady-state  measurements. 
This  ultrafast  short-wavelength  transient  diffraction  methodology  represents  a  newly-developed 
experimental  capability  that  allows  us  to  map  out  nanoscale  thermal  transport  over  a  range  of 
nanostructure  sizes. 


SUMMARY 

To  close,  we  have  made  significant  progresses  over  the  past  3  years  for  the  development  of 
nanostructured  materials  with  tunable  thermal  properties.  We  have  developed  a  novel 
computational  framework  for  the  analysis  and  design  of  nano-structured  materials  with  tunable 
thermal  properties.  This  novel  multiscale  modeling  approach  is  based  on  the  Boltzmann  transport 
equation  (kinetic  theory).  We  have  developed  nonequilibium  Green’s  function  approach  for 
studying  phonon  transmission  across  interfaces  of  dissimilar  materials.  Mechanical  stress  has 
been  explored  as  a  mechanism  for  tuning  thermal  conductivity  of  materials.  We  have  also 
developed  optical  pump  soft  X-ray  probe  based  nano-scale  metrology  for  measuring  nanoscale 
thermal  transport.  Our  work  made  the  first  time  observation  and  quantitative  measurement  of 


non-Fourier  heat  dissipation  from  a  nanoscale  hot  spot  to  the  surroundings.  We  have  also 
identified  future  research  directions  for  the  development  of  nanostructured  materials  with  tunable 
thermal  properties. 


Figure  6:  (a)  Schematic  illustrating  the  difference  between  diffusive  and  quasi-ballistic  thermal 
transport  across  a  nano-interface,  (b)  Sample  geometry  showing  laser  heating  and  EUV 
detection  of  heat  flow  from  nanostructures  to  the  bulk  substrate  underneath,  (c)  Measured 
thermal  resistivity  from  Ni  nanostructures  to  fused  silica  (blue)  or  sapphire  (red)  substrates.  The 
blue  and  red  dashed  curves  show  model  predictions  for  the  ballistic  resistivity  correction, 
assuming  phonon  mean  free  paths  of  2nm  (FS)  and  120nm  (sapphire).  For  linewidths<  200nm 
on  sapphire,  heat  flow  into  the  bulk  is  significantly  slower  than  predicted  by  the  Fourier  Law. 
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